Properties of a fractional derivative Schrodinger type wave equation and a new 
interpretation of the charmonium spectrum 



Richard Herrmann 

GigaHedron, Farnweg 71, D-63225 Langen, Germany 

(Dated: February 7, 2008) 

Based on the Caputo fractional derivative the classical, non relativistic Hamiltonian is quantized 
leading to a fractional Schrodinger type wave equation. The free particle solutions are localized 
in space. Solutions for the infinite potential well and the radial symmetric ground state solution 
are presented. It is shown, that the behaviour of these functions may be reproduced by an ordi- 
nary Schrodinger equation with an additional potential, which is of the form ~ |a;| for a < 1, 
corresponding to the confinement potential, which is introduced phenomenologically to the stan- 
dard models for a non relativistic interpretation of quarkonium-spectra. The ordinary Schrodinger 
equation is triple factorized and yields a fractional wave equation with internal SU{3) symmetry. 
The twofold iterated version of this wave equation shows a direct analogy to the derived frac- 
tional Schrodinger equation. The angular momentum eigenvalues are calculated algebraically. The 
resulting mass formula is applied to the charmonium spectrum and reproduces the experimental 
masses with an accuracy better than 0.1%. Extending the standard charmonium spectrum, three 
additional particles are predicted and associated with Ec(2455) and y(4260) observed recently and 
one X(4965), not yet observed. The root mean square radius for Ec(2455) is calculated to be 
<r >~ 0.3[fm]. The derived results indicate, that a fractional wave equation may be an appropriate 
tool for a description of quark-like particles. 

PACS numbers: 12.39, 12.40, 14.65, 13.66, 11.10, 11.30, 03.65 



I. INTRODUCTION 

Since Newton^] and Leibniz 2] introduced the con- 
cept of infinitesimal calculus, differentiating a function 
f{x^, a;") with respect to the variable is a standard 
technique applied in all branches of physics. The deriva- 
tive operator di, 



d 



(1) 



transforms like a vector, its contraction yields the 
Laplace-operator, using Einstein's sum convention 



(2) 



which is a second order derivative operator, the essential 
contribution to establish a wave equation, which is the 
starting point to describe several kinds of wave phenom- 
ena. 

Until now, in high energy physics the derivative op- 
erator has only been used in integer steps. We want to 
extend the idea of differentation to arbitrary, not neces- 
sarily integer steps. A natural generalization is to search 
for an operator Di by setting 



(3) 



where m, n are integers. Formally, this is solved by ex- 
tracting the m-th root 



D, = a 



71 /m 



(4) 



or, even more general, we will introduce a fractional 
derivative operator by 



a £ 



(5) 



with the fractional derivative coefficient a being a posi- 
tive, real number. 

The concept of fractional calculus has stimulated 
mathematicians since the days of Leibniz Q-Q- In 
physics, early attempts in the field of applications was 
studies on non-local dynamics, e.g^. anomalous diffusion 
or fractional Brownian motion [gjjH- 

During the last decade, remarkable progress has been 
made in the theory of fractional wave equations, 
Raspini,8i,;9| has derived a fractional (a = 2/3) Dirac 
equation. Baleanu[lo| has studied the Euler-Lagrange 
equations for classical fields and gave the explicit form 
of a fractional Klein-Gordon-equation and a fractional 
Dirac-equation, conformal with Raspini's. 

Both studies were based on the use of the Riemann- 
Liouville (RL) fractional derivative, which is used by 
many authors working on the field of fractional deriva- 
tives. 

For practical purposes, the main deficiency of the RL 
fractional derivative is the fact, that the derivative of a 
constant function does not vanish. Maybe this is one 
reason for the fact, that until now, there exists not a 
single ap plic ation in the area of high energy physics. 

Laskinjig has proposed a hermitean fractional 
Schrodinger equation, based on Feynman's path integral 
approach. His applications are based on the semi classi- 
cal Bohr-Sommerfeld quantization condition only. 

We will use a different approach. We will apply the 
concept of fractional derivative to derive a fractional 
Schrodinger type wave equation by a quantization of the 
classical non relativistic Hamiltonian. We will collect ar- 
guments and results which indicate, that this equation 
is an alternative tool for a appropriate description of the 
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charmonium spectrum, which is normally described by a 
phenomenological potential. 

In the following sections, we will explicitely derive ex- 
act solutions for the free particle and for particles in an 
infinite well potential. We will prove, that these solu- 
tions show a behaviour, which may be reproduced by an 
ordinary Schrodinger equation with an additional linear 
potential term for a < 1, indicating that a fractional 
wave equation and the confinement problem are strongly 
related. 

We will derive a fractional multi-component wave 
equation via threefold factorization of the ordinary non 
relativistic Schodinger equation, which contains an inter- 
nal SU{3) symmetry. 

We will then study an analytical mass formula in terms 
of angular momentum multiplets which will reproduce 
the experimental masses of the charmonium spectrum 
within an error of better than 0.1%. 

We will extend the standard charmonium spectrum 
and predict new, additional particles. 

Finally, we will give a reasonable estimate for the root 
mean square radius of 11° (2455). 



II. FRACTIONAL DERIVATIVE 

Let q = [a] be the integer part of a and / a function 
of n variables x^, i — 1, ...,n with > 0. To derive a 
specific representation of the fractional derivative opera- 
tor Di, defined by 10) we start with the Cauchy integral 
extended to fractional order 



r(a) Jo 



(6) 



a formal split of the partial differential operator into a 
fractional integral and an integer differential part 



a _ Qa-(g+l)Qg+l 



9: 



(7) 



leads to the definition of the Caputo 17] fractional differ- 
ential operator ^Z?, which is the form we will use. 



1 du d'i+^ 



(8) 



T{q+l- a) Jo (x* - u)""? 



For a constant function this fractional derivative van- 
ishes: 



^D{x) const = 
For a function of the type 
fix) = x"" 



n e 



(9) 



(10) 



the fractional derivative is: 

^ ' r(l + (n - l)a) ^ ' 

For k > 0, a; > we are then able to define Caputo- Taylor 
series of the form 



f{kx) = ^ a^ikxY' 



(12) 



The corresponding fractional derivatives are given by: 

T{l + {n + l)a) 



(1 + no) 



-(fcx)"" (13) 



Since we intend to use x and the fractional derivative 
D{x) on R, the next step is an extension of our defini- 
tion of X and '^Dix) to negative reals ]R_. 

We propose the following mappings for x — > x(a;) and 
tD{x) D{x): 



X{x) = sign(x) 

D{x) = sign{x)tD{\x\) 



(14) 
(15) 



Besides a unique mapping from to M.^ and M_ to 
the behaviour under parity transformations 11 is well 
defined: 



IlD{x) = -D{x) 



(16) 
(17) 



With these definitions we are able to define series / on 



/(x(M) = E«»^"(^^) 



(18) 



n=0 



with a well defined derivative D 



Dfixikx)) - sign(fc)|fcr E a»+i ^^lT/"lA"^ x"(M 



71=0 



r(l + na) 



(19) 

To construct a Hilbert space on functions f{x{kx)) we 
first define the integral operator 

du" = sign(x)+/"(|a;|") (20) 

with the fractional scalar product 

< f\g / du" rixiku)) gixik'u)) (21) 



An expectation value <0> of an operator O may con- 
sequently be defined with 

<f\0\g>^ [ duTixiku))dg{x{k'uj) (22) 
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to be 



<0>= 



<f\0\g> 



(23) 



<f\9> 

The space coordinates and corresponding derivatives, de- 
fined by ((njl and lfTB|) , are the basic input for our deriva- 
tion of a fractional non relativistic Schrodinger type wave 
equation, the corresponding angular momentum opera- 
tors and caculation of expectation values. 



Thus, a time dependent Schrodinger type equation for 
fractional derivative operators results 



i/"* = (--mc' ( — 
^2 \mc 



2a 



D'D, (34) 

y(x\...,x\...,x3^))* = ihdt-^ 



For a = 1 this reduces to the classical Schrodinger equa- 
tion. 



III. QUANTIZATION OF THE CLASSICAL 
HAMILTONIAN AND FREE PARTICLE 
SOLUTIONS 

By use of the definitions H14(l , 1)15(1 for the space coor- 
dinate and for the fractional derivative, we are able to 
quantize the Hamiltonian of a classical non relativistic 
particle and solve the corresponding Schrodinger equa- 
tion. 

We define the following set of conjugated operators on 
an euclidean space for N particles in space coordinate 
representation: 



= {Po,P^}^{ifl^u-^ 



mcDi} (24) 



{ihdt, — ] mc sign(a;*) +i:'(|a;*|)} (25) 
mc / 



{Xo, Xi} — {t, 



(l-a) 



■} (26) 



(l-a) 



1 



mc J r{a + 1) 



r(a + 1) 
sign(a;,)|a;.r} (27) 



1,...,3N 



These operators satisfy the following commutator rela- 
tions on a frmction set {a;""}: 



Xi , Xj 



X,,P, 



-ih5. 



1 



T{l + na) 



(28) 
(29) 
(30) 



r(l + (n+l)a)) 



(31) 



T{l + {n-l)a) r{l + na) 

= —ihSij c{n, a) 

With these operators, the classical, non relativistic 
Hamilton function He, which depends on the classical 
momenta and coordinates {pi,a;*} 

3N 2 

H,^y2^ + V{x\...,x\...,x''') (32) 
— ' Zm 

i=l 

is quantized. This yields the Hamiltonian 

1 / h \ 

H°' ^ mc^ ( — ) + V(X\ X\ X^^) 

2 \ mc J 

(33) 



A. Properties of the momentum operator 

We extend the standard series expansion of the expo- 
nential function to the fractional case 



exp(a,x(M) = 2^ — 



n=0 



^(1 + an) 



(35) 



^ I7,™|2an 



71- 

= (B„(x(M) 



\kx\ 



(2n+l)Q 



r(l + 2an 

n=0 ^ 



^-r(l + (2n+l)a) 



(36) 

where is the Mittag-Leffler functionfisf . The func- 
tions t/j — exp(Q;, —ix{kx)) are eigenfunctions of the mo- 
mentum operator with the real eigenvalues 



P^ = mc sign(fc) \k\" tp 

mc ' 



(37) 



The Leibniz product rule, which plays an important 
role for the standard derivative, is not valid any more 
for the fractional derivative. Instead, with an arbitrary 
additional function i?, we can write: 

D{fg)^{Df)g + f{Dg) + R (38) 

For the momentum operator it follows 



du'^riPg) = fg- J du'^iPHg- J du"R 

= + [ du'^{Pf)*g~ I du^'R (39) 



Consequently, neither P nor the fractional Schrodinger 
operator with H°' in H34|l are hermitean operators. 

A direct consequence is the non orthogonality of the 
calculated eigenfunctions. 

In general, hermitean operators are preferred, since 
their eigenvalues and expectation values are always real. 
Nevertheless, eigenvalues for momentum, energy and an- 
gular momentum as well as expectation values derived 
with the proposed fractional Schrodinger equation H34() 
turn out to be real, as will be demonstrated in the fol- 
lowing sections. 

Anyhow, we doubt, that a fractional operator should 
always be hermitean. A typical example was the expecta- 
tion value of the root mean square radius of a free quark. 
Indeed, any real value would be a contradiction to the 
experiment. 
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FIG. 1: Free particle solutions for the fractional derivative 
operator Schrodinger type equation, for different values of a. 
Upper part of figure are cos(Q,a;), lower part are sin(Q:,3;) 
with a = 0.85, 0.90, 0.95 (solid thin line), . a = 1 (solid thick 
line) and a = f .05, 1.10 (dashed thin lines) each. For a = 1 
solutions reduce to the standard cos(a;) and sin(a;) functions. 
For a < 1 these functions are increasingly located at a; = 0. 
For a > 1 the amplitudes increase with increasing x. 



FIG. 2: Zeroes of eigenfunctions cos(q:, ^x) (solid lines) and 
sin(a, ^x) (dashed lines) of the free particle eigenfunctions. 
For a < 1/2 there are no zeroes. For a — 1 the zeroes are 
given by fe„ = n € N. For a > 1 an infinite number of zeroes 
exists. 



following relations are valid: 



B. Free particle solutions 

We will now present the free particle solutions for the 
fractional Schrodinger type equation (|34|) . Wc can do 
this, since for V{X'^) = the commutator [P^,H°'] van- 
ishes and consequently, energy and momentum are con- 
served. 

Let us first consider the one dimensional case. We 
extend the definition for the standard series expansion 
for the sine and cosine function 



cos(q;,x) — 



(-l)"l^l 

T{2na + 1) 
(S2a(-x'(fc.x)) 



(40) 
(41) 

sin(a,x) = sign(:r) ^ i^^:^^l-L— ^ (42) 



^^^r((2n+l)a + l) 

X{kx) €2a,l+a{-x'^(kx)) 



(43) 



where and G;^,^ are Mittag-Leffler 18] and generalized 
Mittag-Leffler functions 19]. With these definitions, the 



Dsm{a,kx) = sign(fc)|fc|" cos(a, fcx) (44) 
Dcos{a,kx) = — sign(fc)|fc|" sin(a, /cx) (45) 

It follows from these relations, that the above functions 
(|40|l are the eigenfunctions of the free Schrodinger type 
equation H34|l in one dimension. In the stationary case 
we get the energy relation 



1 



2 \ mc J 



(46) 



This result may easily be extended to the n-dimensional 
case. 

In figure ^ the functions sin(Q;, x) and cos(q;, x) are 
plotted for different values of a. While for a = 1, these 
functions reduce to the known cos(x) and sin(a;), which 
are spread over the whole x-region, for a < 1 these func- 
tions become more and more located at a; = and oscilla- 
tions are damped, a behaviour, which we know e.g. from 
the Airy-functions. For a > 1 the functions amplitude 
increases for increasing x. For a < 1 these functions are 
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normalizable on R: There exists an upper bound M with 

/oo 
dx" cos(a, fca;) cos(a, fc'x) < M (47) 
-OO 

/OO 
dx°' s\n{a,kx)sm{a,k'x) < M (48) 
-OO 

dx°' cos(a, kx) sin(a, k'x) — (49) 



For a > 1 this integral is not bound any more, instead 
a box-normaUzation with a box size L much larger than 
the dimensions of the system considered is proposed. 



C. Particle in an infinite potential well 

Now we will give the exact eigenfunctions and eigen- 
values for a particle confined in an infinite potential well. 
We first investigate the one dimensional case. Therefore 
we define the potential V{x) 



V{x) 



\x\ < a 
\x\ > a 



(50) 



The corresponding boundary condition for the eigenfunc- 
tions ^p{x) is: 

i/;(-a) = ^(a) = (51) 

In figure |3 the zeroes for the free particle solutions 
cos(a, ^x) and sin(a, ^x) are plotted. For a < 1/2 there 
are no zeroes. For the interval 1/2 < a < 1 exists only 
a finite set of zeroes. For a > 1 an infinite number of 
zeroes exists. 

Let be a zero of the free particle solutions. The 
eigenfunctions of the infinite potential well potential are 
then given by: 

(52) 



4l\x) = cos{a,x{kL-)) 
V'Lli(a^) = sin(a,x(fc2«+i7)) 



(53) 



where the sign indicates the parity of the states. The 
normalization condition is 



(54) 



dx"^*^{x)^pn{x) = 1 

The energy is then given by 

2 \mc / a 
The extension to the N-dimensional case is then 

N 

and for the energy 

i5n.„« = ^( — ) mc^El — I'" 
2 V mc / ^ — ' " ■ 



i=l 



\ f fi \ o X — A . k^ 



(55) 



(56) 



(57) 



a=0.9 
n=6 A 






n=5Vy 




n=4\ / 
n=3 / 









a=1.1 



n=6 




n=5 








FIG. 3: The six lowest eigenfunctions for the one dimensional 
infinite square well potential, with a — 0.9 on the left and 
a = 1.1 on the right side. Below the corresponding potential 
for equivalent solutions of the ordinary (q = 1) Schrodinger 
equation, according to l|71|l is plotted. For a < 1 this potential 
contains a strong linear admixture, for a > 1 to potential 
graph behaves like ~ 1 — a;^ 
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D. Radial solutions 

In the case of fractional derivative operators there ex- 
ists no general theory of covariant coordinate transfor- 
mations until now. 

We intend to perform a coordinate transformation 
from carthesian to hypersphcrical coordinates in 

f{xi,X2,...,XN) = /(r, 01,02, ••■,(?!>7V-i) (58) 

The invariant line element in the case a = 1 

ds^ = g,jdx'dx^ i,j = 1, ...,N (59) 

for arbitrary fractional derivative coefficient a may be 
generalized to 

ds^" = g'^jdx'^dx^'^ i,j = l,...,N (60) 

Consequently a natural definition of the radial coordinate 
is given by 



N 



(61) 



We assume the spherical ground state to be independent 
of the angular variables, square integrable and of positive 
parity. Therefore an appropriate ansatz is 

oo 

g{N, a, kr) = ^(-l)"a„(^, a)(|fc|r)2™ (62) 

n=0 



or in carthesian coordinates 



N 



g{N, a, fcxi, ...,kxN) = ^(-l)"a„(iV, a)(^ jfexj^")" 

n=0 1=1 

(63) 

where the coefficients a„ depend on the explicit form of 
the potential. 

For a free particle, a solution on is given with the 
abbreviation 

7?, = r(l + 2aj)/r(l + 2a{3 - 1)) (64) 

by the recurrency relation 

flo = 1 

aj = aj^il {{N + j = 1,2,... (65) 

An infinite spherical well is described by the potential 



V{r) 




(66) 



The corresponding boundary condition for the ground 
state wave function g{N, a, r) is: 



Let fcgpjj be the first zero of the free particle ground state 
wave function, the ground state wave function for the 
spherical infinite well potential is given by 

g{N,a,k%^r/ro) (68) 

and the ground state energy is then given by 



2 \ mc tq 



(69) 



E. 



Remarks on equivalent solutions for the 
ordinary Schrodinger equation 



We have shown, that the free particle solutions and 
the solutions for the infinite potential square well of the 
fractional Schrodinger equation for a < 1 are localized at 
the origin and for a > 1 are localized at the boundaries 
of a given region respectively. 

We want to deduce a similar behaviour of these func- 
tions in terms of the ordinary [a = 1) Schrodinger equa- 
tion. Let us assume, the eigenf unctions of the fractional 
Schrodinger equation may be equivalently interpreted as 
solutions of the ordinary (a = 1) Schrodinger equation 
with an additional potential V . 

In order to derive the explicit form of this potential, we 
use the following relation between the eigenfunctions 
temperature T and a given potential V ^ which is derived 
within the framework of thermodynamics and statistical 
quantum mechanics [25l|: 

Let '^n and En be the eigenfunctions and energy eigen- 
values of the ordinary Schrodinger equation with a given 
potential V: 



(70) 



— A + y ) ^n=En^n 



As long as the temperature T is large compared to the 
average level spacing, the relation [N is the normalization 
constant) 



-^exp(-y/T) ^ 



(71) 



g(iV,a,ro) = 



(67) 



T » Ai?„ — En+l — En 

is valid. 

Since eigenfunctions and eigenvalues for the fractional 
free particle solutions are known, we therefore are able 
to deduce the explicit form of such a potential. 

In figureOlthe graph of V/T is plotted for a — 0.9 and 
a = 1.1 respectively. For a = 0.9 the potential contains 
a dominant linear term V ^ \x\ , while for a = 1.1 a 
behaviour like V ^ 1 — may be deduced. Of course, 
for a = 1 a potential V = const results. 



The behaviour of the fractional eigenfunctions may al- 
ternatively be interpreted, assuming that a is a mea- 
sure of charge for a particle moving in an homogenously 
charged background. 

For a > 1, we could assume a charged particle, with 
charge e.g. a moving in a homogeneous background of 
charged particles with charge 1 — a. Since for a homoge- 
neously charged sphere and approximately for a charged 
box too, the potential inside the box is ~ ^ ''^) 

with Z — 1 — a. This obviously would explain a repulsive 
force. The energetically favoured positions for a particle 
with charge a indeed were the boundaries of the box. 

For a = 1 the background was neutral and therefore 
no additional interaction with a particle was present. 

For a < 1 this simple model could explain an attrac- 
tion, but not the linearity of the potential. 

Therefore we obtain the remarkable result, that in 
the case a < 1, the free particle solutions of the frac- 
tional Schrodinger equation show a behaviour, which is 
equivalent to the behaviour of solutions of the ordinary 
Schrodinger equation with an additional linear potential 
term. In other words, the solutions of a free fractional 
wave equation with a < 1 automatically show confine- 
ment, a property, which was first observed for quarks. 

In order to obtain more properties of the fractional 
derivative operator Schrodinger equation, we will now 
calculate the eigenvalues of the angular momentum op- 
erator. 



IV. CLASSIFICATION OF ANGULAR 
MOMENTUM EIGENSTATES 

We define the generators of infinitesimal rotations in 
the i,j-plane = 1, ...,3A'^), with N being the number 
of particles): 



V{a + l) 



r(a + 1 



■A 



(72) 



We will derive the angular momentum eigenvalues alge- 
braically. Thus it is necessary, to apply the commutator 
relation [^i,Pj] (see repeatedly. In figure^ this 

commutator c(n, a) is plotted. It shows a smooth depen- 
dence on n, which we have to eliminate. 

Ignoring both, the n and a dependence we set as a 
lowest order approximation: 



Co 



c(0,l) = l 



(73) 



A more precise statement for c(n, a) can be deduced from 
the following consideration: Since we will concentrate on 
the lowest enery levels only , the approximation 



Ci{a) — c(0, a) = 1 — 



1 



r(i -a)r(i 



a] 



(74) 



is valid. Of course, this overestimates the commutator 
for higher values of n. Therefore a more sophisticated 




n 



FIG. 4; Commutator [Xi,Pi\ from H^O^ in units of h on the 
function set {a;""} for different values of a. The dependence 
on n is a direct consequence of the fact, that the Leibniz 
product rule is not valid any more for fractional derivatives. 
For a = 0.6 the positions of the successive approximations 
on c(n, a) according to II73II . II74II and 17511 are given as circle, 
triangle and squares. 



treatment fixes c{n, a) for a given j to be: 

F(l + (j + l)a) F(l+ja) 



C2(i,a) = 



F(l + ja)T{l + a) r(l + (j - l)a)F(l + a) 

(75) 

We will consider these three cases, which allows an es- 
timate on the validity of the results. Therefore, as long 
as the commutator does not depend on n, [Lij, H°'] van- 
ishes and angular momentum is conserved. Commutator 
relations for Ly are isomorph to an extended fractional 
S'0"(37V) algebra: 



_(76) 

Consequently, we can proceed in a standard way [^^ , by 
defining the Casimir operators 



.,3iV (77) 



which indeed fulfill the relations [Ag^jLy] = and suc- 
cessively [A^,A^,] = 0. Consequently the successive 
group chain 

SO°'{iN) D 50"(3iV - 1) D . . . D S'0"(3) D SO°'{2) 

(78) 

is established. The explicit form of the Casimir operators 
is given by 

hi ^ +X'X,P^Pj - ihc{a) {SI - l)X^Pj - X'X^P.Pj 

(79) 

We introduce a generalization of the homogeneous Euler 
operator for fractional derivative operators 



J,^(a)=x(a:*)/F(l + a)A 



(80) 
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TABLE L Eigenvalues in units of h for the angular momentum states of a single particle, n is a counter for the eigenvalues 
of the Euler operator, eigenvalues for L2{a) and Jc{a) are given for a — 1, a = 2/3, Q = 0.68 and a. = 0.65. J|(a) are listed 
with c according to 1731 1. 1741 1 and 17511 . 





iz(2/3) 


Lz(0.68) 


^o'(l) 


^'(2/3) 


J|(0.68) 


J? (0.65) 


J|(j, 0.65) 


























1 


1 


1 


2 


2 


2 


1.604 767 


1.478 157 


2 


1.460 998 


1.478 157 


6 


3.595 515 


3.663 108 


3.078 892 


2.800 590 


3 


1.860 735 


1.894 649 


12 


5.323 069 


5.484 346 


4.735 519 


4.305 776 


4 


2.222 222 


2.272 597 


20 


7.160 493 


7.437 298 


6.539 094 


5.961 779 


5 


2.556 747 


2.623 332 


30 


9.093 704 


9.505 205 


8.468 379 


7.747 796 


6 


2.870 848 


2.953 417 


42 


11.112 618 


11.676 094 


10.508 808 


9.649 033 



With the generahzed Euler operator the Casimir- 
operators are: 

Kl = +X^X,PP, + ft2 (k - 2)J^ + J^J^) (81) 

Now we define a Hilbert space ^q, of all homogeneous 
functions /, which satisfy the Laplace equation D^Dif — 
and are normalized in the interval [—1,1]: 

i3a = {/ : /(Af) - A""/(a;); &DJ = 0} neN 

(82) 

This is the quantization condition. It guarantees, that 
solutions are regular at the origin. 

On this Hilbert space, the generalized Euler operator 
Jg (a) is diagonal and has the eigenvalues 



lk{a,n) = 







r(na + 1) 



r(a + l)r((n-l)a + l) 



for rt = 
for n — 1,2, 3, ... 

(83) 

This is the main result of our derivation. 

We want to emphasize, that these eigenvalues are dif- 
ferent from the degree of homogenity in the general case 
a 7^ 1, or, in other words: only in the case of a = 1 the 
homogenity degree n of the polynoms considered coin- 
cides with the eigenvalues of Jg (a = I,")- 

Once the eigenvalues of the generalized Euler opera- 
tor are known, the eigenvalues of the Casimir-operators 



A2, A? are known, too: 



A2/ 



hl2{a,n)f 
f?lk{a,n){lk{a,n) 



c(a) (fc-2))/ 



with 



Ih. > lu 



fe-i 



> 



>l ±h l> 



(84) 
(85) 



(86) 



For the case of only one particle [N — 1), we can intro- 
duce the quantum numbers j and m, which now denote 
the j-th or m-th eigenvalue of the Euler operator. The 
eigenf unctions are fully determined by these two quan- 
tum numbers / = \jm> 

With the definitions = L12 and 



^23 



it follows 

Lz \ jm> = hl2{a,m) \jm> (87) 
m = 0,±l,±2,...,±j 

J2|jm> = h%{a,j){h{a.j) + c{a))\jm> (88) 
j = 0,-^1,-1-2,... 

Please note the fact that Lz remains unchanged for any 
choice of constant c(a), only changes. 

In tableQlthe first seven eigenvalues of Lz and for a 
single particle are listed for a = 1, a = 2/3,0.68 and 
a = 0.65 and different approximations for c(a). For 
a ^ 1 the eigenvalues of the generalized Euler operator 
are not equally spaced any more. For a < 1 the stepsize 
is strongly reduced. Since the generalized Euler opera- 
tor eigenvalues contribute quadratically into the defini- 
tion J^, the energy of higher total angular momenta is 
reduced increasingly. 

We have derived the full spectrum of the angular mo- 
mentum operator for the fractional derivative operator 
Schrodinger type wave equation by use of standard alge- 
braic methods. 

We will get additional information about the proper- 
ties of this wave equation, if we consider its factorized 
pendant. We present some results in the next section. 



V. RESULTS FOR THE FACTORIZATION OF A 
NON RELATIVISTIC SECOND ORDER 
DIFFERENTIAL EQUATION 

Linearization of a relativistic second order wave equa- 
tion was first considered by Diracj^J- Starting with 
the relativistic Klein-Gordon equation his derived Dirac 
equation gave a correct description of the spin and the 
magnetic moment of the electron. 

The concept of linearization is important, since it pro- 
vides a well defined mechanism to add an additional 
SU(2) symmetry to a given set of symmetry properties 
of a second order wave equation. 

Since linearization may be interpreted as a special case 
of factorization, namely to 2 factors, a natural general- 
ization is a factorization to n factors. 
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In 2000, Raspinifs'l proposed a Dirac-like equation with 
fractional derivatives of order 2/3 and found the corre- 
sponding matrix algebra to be related to generalized Clif- 
ford algebras; in 2002 Zavada generalized Dirac's ap- 
proach, and found, that rclativistic covariant equations 
generated by taking the n-th root of the Klein-Gordon 
or d'Alembert operator ( D^/") are fractional wave equa- 
tions with an additional SU(n) symmetry. 

These results indicate, that fractional order wave equa- 
tions may be appropriate candidates for a description of 
particles, which own a SU(n) symmetry. The case n = 3, 
which corresponds to a triple factorization is therefore 
important for a description of particles with a SU(3) sym- 
metry. 

Whether or not a factorization of non rclativistic wave 
equations leads to similar results, has not been examined 
yet. 

In 1967, Levv-Leblond (2^ has linearized the non rcla- 
tivistic Schrodinger equation and obtained a linear wave 
equation with an additional SU(2) symmetry, but until 
now his approach has not been extended to higher frac- 
tional order. 

In order to obtain additional information on the inher- 
ent symmetries of the fractional Schrodinger equation, 
which we proposed in (|34|l , we therefore derive in the fol- 
lowing section the explicit form of a fractional operator, 
which evolves from a triple factorization of the ordinary 
Schrodinger equation: 



A. Triple factorization of the non relativistic 
Schrodinger equation 

We intend to derive a fractional operator R, which, 
iterated 3 times, conforms with the ordinary, non rela- 
tivistic Schrodinger operator: 



RR'R" 



2m 



-A - ihdf 1 



(89) 



where In is the n x n unit matrix. 
We use the following ansatz: 



R 
R' 
R" 



= aAdI' + bB'd^' + cC 



aA"d^* + bB"'d^^ + cC" 



(90) 
(91) 
(92) 



with matrices ^, A', A", B, B', B", C, C", C", fractional 
derivative coefScients for time and space derivative at , 
and scalar factors a, 6, c, which will be determined in the 
following. According to Zavada T^, we define a triad of 
unitary, traceless 3x3 Pauli type matrices, which span 
a subspace of SU (3) with 



exp(^ k) 



an explicit representation is 

xi 0\ 

(94) 
(95) 
(96) 

These matrices obey an extended Clifford algebra 

a' a' a'' = 6 S'^'' i, j, fc = 1, 2, 3 (97) 

all Permutations 

Let (g) denote the outer product of any two matrices. In 
order to describe a single particle with the coordinates 
{t, x,y, z}, we define the following 4 7 matrices, with 
dimension 9x9: 











n 


n 


X2 


3^3 


n 


n 





X2 











Xl 










Xi 






















X3 



7 



1.3 «) cr'' 



Y ^ a'(E)a^ i = 1,2,3 



(98) 
(99) 



Now we are able to specify the above introduced matrices: 

(100) 

:(70-X2l9) 



A = -L(70-xil9) 



A' 



1 

1 



(101) 



A" ^ _(70_x3l9) (102) 

B' = 7* (103) 

B'' = 7* (104) 

B"' = f (105) 

/I 0\ 

C = a;il3® (106) 

\o oy 
/O o\ 

C = a;2 13 <^ 1 (107) 

\o oy 

/O o\ 

C" = X3I3® \ o\ (108) 
\q l) 

with these specifications RR'R" yields 

RR'R" = (^a^cdf'" +b'^Yl ^^^""^ ^9 (1^9) 

A term by term comparison with the nonrelativistic 
Schrodinger operator determines the fractional derivative 
coefficients: 



1,2,3 



(93) 



at = 1/2 
a, = 2/3 



(110) 
(111) 
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and the scalar factors: 



= — 



1/6 



2ni 



(112) 

(113) 
(114) 



Finally, according to (|14|l and we extend the deriva- 
tive operator on M via: 



da 



sign{x^)tD{\x^\) 



(115) 



Thus the fractional operators R, R', R" are completely 
determined. As a remarkable fact we note the different 
fractional derivative coefficients for the fractional time 
and space derivative. 

We therefore have proven, that the resulting SU(3) 
symmetry is neither a consequence of a relativistic treat- 
ment nor is it a specific property of linearization (which 
means, specific to first order derivatives of time and space 
respectively) . 

It is a consequence of the triple factorization solely. 

We want to emphasize, that the factorization not only 
determines the symmetry group of the 7-matrices used, 
but also determines the dynamics of the system, forcing 
e.g. for a SU(3) symmetry fractional time {at = 1/2) 
and space derivatives(ai = 2/3). 

Thus, the result of factorization shows two important 
properties: It yields an additional SU(n) symmetry and 
simultaneously the corresponding dynamics. Applied to 
QCD, this seems more consistent than the standard con- 
cept based on Yang-Mills field theories, where an arbi- 
trary non abelian symmetry group gauge field, which first 
has to be deduced from experimental data to be a SU(3) 
symmetry is coupled to a symmetry independent dynam- 
ical (e. g. Dirac) field, neglecting a possible influence of 
the symmetry onto the dynamics. 

In that sense, the fractional operator R, defined in (|90|l . 
would be an alternative starting point for a pure, non rel- 
ativistic QCD, since it contains a consistent description 
of both, symmetry and dynamics of a pure SU(3) sym- 
metry, without any additional SU(2) admixture. 

Finally, besides the fractional wave equation operator 
R and the triple iterated RR'R", which corresponds to 
the ordinary Schrodinger operator, an additional type of 
wave equation, the twofold iterated RR' emerges, which 
reads: 

cRR' = a^cAA'dt + b^cB'B' df'^d^/^ + additional terms 

(116) 

or, inserting the factors: 

5(2) = _^MA'a,_i(A)4/3^^2( 1)1/3 5.5,52/3.2/3 

2 mc 2 ' 



-additional terms 



117) 



If) 
(/) 

CO 



i500 








\[/ (4410); <40> ^— 
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prediction 2 Y(4259): <33> 


m=3. 




i|/(4160);<32> 


m=2 




»|/ (4040); <31> ^— 


m=1 




\|/ (J / iV), <j\j> ^^^^ 


m~0 




X (3bbb), <iL^> 


m=2 


3500 


X(3511);<21> 


m=1- 




X(3415);<20> 


m=0 




j;v(3097);<11> 


m=1 


3000 


il„(2980);<10> 


m=0- 




Standard spectrum 


2500 








(2451); <oo> prediction 1 


m=0 




j=0 j=1 j=2 j=3 j=4 





Obviously S"^^) and the fractional Schrodinger equation 
(|34|l . we derived in section ITTll are closely related. 



FIG. 5: Charmonium spectrum. All observed particles are 
given with their name, experimental mass from|2q| and the 
proposed SO'^i?)) conforming quantum numbers j and m, 
which are the j-th and m-th eigenvalue of the generalized Eu- 
ler operator. Predicted particles are: in the lower part of 
the spectrum denoted as prediction 1 < 00 > associated with 
E1?(2455) and in the upper part prediction 2 <33> associated 
with Y(4260). 



i?" seems nothing else but a scalar version of 5(2) for 
the special case a = 2/3. Therefore, examination of the 
properties of the scalar fractional Schrodinger equation 
(|34|l with a = 2/3 should reveal some properties of an 
inherent SU (3) symmetry. 

Summarizing all facts collected, we assume, that the 
fractional derivative operator Schrodinger type wave 
equation with a = 2/3 is an appropriate candidate for 
a non relativistic description of particles with quark-like 
properties. 



VI. INTERPRETATION OF THE 
CHARMONIUM SPECTRUM 

In the previous sections we have introduced the con- 
cept of fractional derivative operators and discussed 
some properties of the resulting non relativistic fractional 
Schrodinger type equation and its factorized pendant. 

We have developed a new theoretical concept, which 
fulfills at least the following three demands: First, avail- 
able experimental data will be reproduced with a reason- 
able accuracy. Second, it will give new insights on un- 
derlying symmetries and properties of the objects under 
consideration. Third, we will make predictions, which 
can be proven by experiment. 
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None of our results, presented so far, does require any 
information from QCD or a similar theory. All our state- 
ments could have been made in the 1930s already, even 
though they would have been highly speculative. To- 
day, we arc in the comfortable position, that there are 
enough experimental data, our predictions can be com- 
pared with. 

A promising candidate is the charmonium 
spectrum p^. 

In the upper part of figure [31 we have displayed all 
experimentally observed charmonium-states with experi- 
mental masses, which are normally compared with results 
from a potential model, which tries to simulate confine- 
ment and attraction by fitting a model potential(26j . [27] . 

We will assume, that this spectrum is a single particle 
spectrum for a particle, whose properties are described 
by the free fractional Schrodinger type eauation H34(l . We 
suppose, the system is rotating in a minimally coupled 
field, which causes a magnetic field Bj. This leads to the 
following Hamiltonian H or mass formula 

H{j,m) |jm>::= [nj'^/h'^ + BjLjh + niQC^) \jm> 

where k, Bj and toqc^ will be adjusted to the experimen- 
tal data. 

The eigenfunctions | jm > are modified spherical har- 
monics with D^Di \jm>— 0, the eigenvalues for and 
Lz are given by l|S7|) . and are hsted in table 1. 

A first remarkable observation is the fact, that states 
with Lz < are missing in the experimental spectrum. 
Only right-handed particles are realized. 

This may be due to the fact, that actually is a 
Casimir-Operator of 50" (2), while Lz is not and there- 
fore the multipletts should more precisely be classified 
according to the quantum number in the general case 
a =/= 1. Consequently in the following we will work with 
Lz>0. 

To check the influence of our approximations (|73|l , (|74() 
and H75I) of c{a) given in H3U|) we will proceed in two steps: 
First, we will consider the case cq = 1. In a second step, 
we will test the influence of the successive approximations 
for c(a) on the accuracy of the proposed mass formula 
with least square flts on the charmonium spectrum. 



A. Interpretation of the charmonium spectrum in 

the case Co = 1 



We first consider the case c{a) = 1. The corresponding 
and Lz values used are listed in tabled as Lz{a) and 
J2(a). 

The first crucial test is the verification of the correct 
value of the non trivial m = 2 quantum number which 
corresponds to the n = 2 eigenvalue of the generalized 
Euler operator. 

For the set of ^-particles (j = 2) (experimental masses 



and errors are taken from (Ig), we obtain: 

X(3556) < 22 > -x(3415) < 20 > 



L,(i = 2,m = 2)e 



X(3511)<21> -x(3415)<20> 
1.478 ±0.007 (119) 



Thus, an a from the experimental spectrum is deduced: 

(120) 



al=^ = 0.680 ±0.006 



Which is remarkably close to the theoretically expected 
ath = 2/3. 

An alternative approach to determine the experimental 
value for a or the m — 2 quantum number repectively, 
follows from the set of ^'-particles {j — 3): 



Lz{j = 3, m = 2)cxp = 



^'(4160)<32> -'5(3770) <30> 
*(4040) < 31 > -^'(3770) < 30 > 
1.44 ±0.09 (121) 



A second experimental a is obtained: 



cxp 



= 0.65 ±0.08 



(122) 



Within experimental errors, both values are identical. 
This observation supports the assumption, that the spec- 
trum may be interpreted using one unique a. 

According to our level scheme, the < 33 > state is 
missing in the standard charmonium spectrum. Using 
0.680 and Lz{a,m) from table the predicted 



mass is: 



^'<33> = 



Lz{a,m ^ 3) 
Lz{a, m = 2) 
4268 ±22 [Me V] 



(^'(4160)<32> -«'(3770)<30>) 

(123) 



In June 2005, the Babar collaboration announced the dis- 
covery of a new charmonium state named Y(4260) [sD]- 
The reported mass of 4259 [MeV] is in excellent agree- 
ment with our mass prediction for \1/ < 33 >. Therefore, 
we associate the predicted particle with y(4260). 

Next we will determine the constants mgc^ and k in 

We choose the two lowest experimental states of the 
standard charmonium spectrum, r]c{2980) < 10 > and 
X(3415) < 20 >. With al^^ = 0.680 and J^iaJ) from 
table n we obtain a set of equations 

moc^ ±2k = 7]c(2980)<10> (124) 
moc^ -h 3.663108k = x(3415)<20> (125) 



which determine 



TOoc^ = 2455 ±3 [MeV] 
K = 262.4 ±0.9[MeV] 



(126) 



Our level scheme predicts a particle with quantum 
numbers < 00 >, which is beyond the scope of charmo- 
nium potential models. According to our mass formula, 
it has a predicted mass of <00>= 2455 ± 3[MeV]. 
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Since this is a low lying state, it should already have 
been observed. Indeed, there exists an appropriate candi- 
date, the E°(2455) < 00 > baryon, with an experimental 
mass of 2452.2[MeV]. This is a charmed baryon with 
quark content (ddc). 

The minimal difference of only 2.8 [Me V] between pre- 
dicted and experimental mass of the E°(2455) < 00 > 
particle indicates, that the assumed fractional SO" (3) 
symmetry is fulfilled exactly. 

Obviously, the fractional SO" (3) multipletts describe 
mesonic and baryonic states of the charm-quark simulta- 
neously. 

Due to its experimentally observed properties, the in- 
ternal structure of the F(4260)<33> particle is subject 
of actual discussion's!^. Besides being a conventional cc 
state, it could alternatively be a tetraquark with con- 
stituents (uucc) or a hybrid charmonium. 

Finally we have only one experimental candidate for 
j = 3 and j = 4 respectively. With parameters H126|l 
using the mass formula (|118|l we obtain the theoretical 
values 

«'(3770)<30>tft = 3894±8[MeV] (127) 
*(4410)<40>(,, 4406±10[MeV] (128) 

For j = 3 the calculated mass (jHTjl differs by 124[Mey] 
from the experimental value. 

On the other hand, the theoretical ^'(4410) < 40 > 
mass ()128ll matches exactly with the experimental value 
within the experimental errors. 

This indicates, that the particles for j — 3, observed 
in experiment, carry an additional property, which re- 
duces the mass by the amount of e.g. a pion. Of course, 
if we add an additional (At)(5j3 term to the proposed 
mass formula, we can shift these levels by the necessary 
amount. 

Summarizing these results, the charmonium spectrum 
reveals an underlying 5*0" (3) symmetry, which agrees 
with the predictions of our theory in the case of a « 2/3. 
The eigenvalues of the generalized Euler operator con- 
form within experimental errors with experimental data. 
Extending the standard charmonium spectrum, two addi- 
tional particles have been predicted and associated with 
E°(2455) and r(4260) observed recently. 



B. Least square fits of the charmonium spectrum 

In the previous section we gave an interpretation of the 
charmonium spectrum for the case c(a) = 1. We found, 
that the spectrum may be described quantitatively, using 
the proposed mass formula for a = 0.680. 

Extending the mass formula (|118|l including a correc- 
tion term for the j = 3 multiplett, we use 

H{j, m) \jm >— kJ'^ /h^ + BjLz/h + niQC^ + (^3 ^ Ar | jm> 

(129) 

to find a fit on the experimental charmonium spectrum. 



To prove, that for cq = 1, a = 0.68 indeed is the ap- 
propriate choice for an interpretation of the charmonium 
spectrum, we minimized errors with respect to a and ob- 
tained a — 0.681. In table HTl the optimum parameter 
sets for a = 2/3 and co,ci,C2 and resulting errors are 
tabulated. 

For ci{a) from (|74|l . a = 0.68 is not the optimum 
choice any more. We therefore minimized errors with 
respect to a, finding a — 0.647 for this case. 

For C2{j,a) from H75|l . we observe a minimal shift in 
a, finding a = 0.649 for this case. This indicates, that 
a more sophisticated treatment of c(a) will only cause 
neglible changes for a and corresponding parameter sets. 

Comparing the optimum parameter sets, the changes 
in the treatment of c(a) are mainly absorbed by the pa- 
rameter K. Parameter At remains remarkably constant. 
This supports interpretation for this parameter to be a 
j — 3 specific quantum number. 

A comparison of experimental with calculated masses, 
based on the optimum parameter sets, is given in table 
mil Mass differences are less than 0.1% and decreasing 
for co,ci,C2. 

With the optimum parameter sets we can predict the 
mass of the < 50 > state to be 

X<50>th= 4965± 10[MeV] (130) 

This state has not been observed in experiments yet. 

We conclude, that the proposed 50" (3) symmetry is 
fulfilled exactly. The values of a = 0.68 and a = 0.65 
resulting from the least square fits are close to , but dif- 
fer significantly from the theoretically expected a = 2/3. 
This indicates, that the inherent SU{3) symmetry is al- 
most fulfilled exactly, with a difference of only 2%. 



C. Size estimate for S° 

Up to now we have treated toqc^ as a simple parameter 
of the proposed mass formula. As a result of our discus- 
sion above, we associate H{0, 0) = toqC^ with the mass 
of If H{j,m) was the solution of the Schrodinger 
equation with a given potential V, then H{0, 0) would be 
interpreted as the zero point energy in this potential plus 
the rest mass of its constituents. 

We intend to estimate the size of E°. Therefore we 
will calculate the expectation value < f > of the radius 
operator f, which is given by 



f = ^X^,+Xl+Xi (131) 

= (—) v(^\ J xr+xl»+xl" (132) 
\7nc J 1 (1 + a) V 

For a first estimate, we choose the infinite square well 
potential H50|) . The energy eigenvalues are given by |S7J). 
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TABLE II: Optimum parameter sets for a fit of the experimental charmonium spectrum with mass formula 11291 in units 
[MeV]. Errors Am are given for the subset m — and the full spectrum Amaii- The first row corresponds to a = 2/3 fixed, the 
following three rows correspond to the three approximations for co,ci and C2 according to 1731 1. 1741 1 and 17511 with a optimized. 
In table UTTl resulting theoretical masses and errors are listed. 



a 


c{a) 






Bi 


B2 




At 




Amaii 


comment 


2/3 


1.00 


2439.33 


274.66 


108.25 


87.00 


263.69 


-129.04 


5.68 


8.98 


a fixed 


0.681 


1.00 


2451.26 


263.83 


117.98 


93.39 


259.16 


-124.48 


1.86 


2.02 


a variation 


0.647 


0.545 


2452.67 


336.16 


119.79 


95.72 


270.19 


-129.00 


1.14 


1.15 


a variation 


0.649 


c(.],a) 


2451.90 


367.41 


116.13 


98.37 


269.46 


-124.39 


0.73 


0.79 


a variation 



TABLE III: Comparison of experimental and calculated masses according to mass formula 112911 with optimized parameter sets 
listed in tabled in units [MeV]. The last row lists predicted theoretical masses for X<50>. 









a = 2/3, Co 


= 1 


Q — 0.68, Co 


= 1 


a = 0.647, ci(a) 


Q = 0.649, C2(i,Q) 


<jm > 


symbol 




mth 


Anith 


mth 


Amth 


mth 


Amth 


mth 


Amth 


<00> 




2452.2 


2439.33 


-12.87 


2451.26 


-0.93 


2452.67 


0.47 


2451.90 


-0.30 


<10> 




2979.6 


2988.66 


9.05 


2978.94 


-0.66 


2977.12 


-2.47 


2980.78 


1.18 


<11> 


J/* 


3096.9 


3096.92 





3096.92 





3096.92 





3096.92 





<20> 


Xo 


3415.2 


3426.89 


11.70 


3417.80 


2.60 


3417.15 


1.95 


3413.65 


-1.54 


<21> 


Xi 


3510.6 


3513.89 


3.30 


3511.19 


0.60 


3512.87 


2.27 


3512.03 


1.43 


<22> 


X2 


3556.3 


3554.00 


-2.26 


3555.85 


-0.40 


3554.68 


-1.58 


3555.26 


1.00 


<30> 




3770 


3772.35 


2.34 


3773.92 


3.92 


3770.17 


0.17 


3770.41 


0.41 


<31> 




4040 


4036.04 


-3.95 


4033.08 


-6.91 


4040.37 


0.37 


4039.87 


-0.13 


<32> 




4160 


4157.60 


-2.39 


4157.02 


-2.98 


4158.39 


-1.61 


4158.30 


-1.70 


<33> 


y(4260) 


4259 


4263.01 


4.01 


4264.98 


5.97 


4260.08 


1.07 


4260.42 


1.42 


<40> 




4415 


4406.07 


-8.92 


4413.80 


-1.20 


4414.36 


-0.64 


4415.22 


0.22 


<50> 


X 




4937.06 




4959.54 




4957.54 




4969.07 





Therefore 



(133) 



2 \ nirC / a 



determines the half boxsize a. The wave function was 
defined in H4U|) . Therefore with the abbreviations 

dV" = dx'^dx^dx'^ (134) 
'^{xi,X2,X3) — cos{kQXi/a)cos{kQX2/a)cos{kQX3/a) 

the expectation value according to (|23|l is 



Setting a = 2/3 and 



mdC 



2452. 2 [MeV] 
300 [MeV] 



(135) 



(136) 



TOcC^ = 1400 [MeV] 
= 1.1648 7r/2 

we derive a = 0.81 [fm] and therefore we obtain the ex- 
pectation value for the radius 



<f(I];;)>cubo= 0.32[fm] 



(137) 



Similarly, we can proceed for the infinite spherical well 
potential (|66|) . 

For the spherical ground state wave function in carthe- 
sian coordinates (|63(l . we obtain 



L^giN = 3,a,xi,X2,x3) 
J"^ g{N = 3,a,xi,X2,X3) ~ 

and therefore g indeed is the ground state |00>. 
The ground state energy bq is given by I69|l 

E° = (2mrf + TOc)c^ + eo(Af = 3,a) 



(138) 
(139) 



(140) 



1 / t \ 2a 1,0 

= (2m, + me)c2 + i — rn,c'\^\'" 

2 \mccj ro 

with fcgpj^ = 3.1652 7r/2 this determines ro = 1.08[fm]. 

With ()135|l we obtain for the expectation value of the 
radius 



<f(E^)>,phcrc=0.33[fm] 



(141) 



Consequently, both potentials lead to similar expectation 
values. 

Since Sj! is not within the scope of standard charmo- 
nium models, there is no direct comparison. 

Nevertheless, there are radii, derived from charmonium 
model calculations p^. reported for |32j]-[3^: 

<f{J/iJj<ll>)> w 0.2[fm] (142) 
<f(xo<20>)> w 0.3[fm] (143) 
<f(4'<30>)> w 0.4[fm] (144) 
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Therefore our results are reasonable compared with 
these calculations. 



VII. CONCLUSION 

Based on the Caputo fractional derivative, we have de- 
fined a fractional derivative operator for arbitrary frac- 
tional order a. A Schrodinger type wave equation, de- 
rived by quantization of the classical non relativistic 
Hamiltonian, generates free particle solutions, which are 
confined to a certain region of space. Therefore confine- 
ment is a natural consequence of the use of a fractional 
wave equation. 

The multiplets of the generalized angular momentum 
operator have been classified acoording to the SO" (3) 
scheme, the spcctnmi of the Casimir-Operators has been 
calculated analytically. 

We have also shown, that for a = 2/3, corresponding 
to a fractional non relativistic Levy-Leblond wave func- 
tion an inherent SU(3) symmetry is apparent. 

Prom a detailed discussion of the charmonium spec- 
trum we conclude, that the spectrum may be understood 
quantitatively within the framework of our theory. Ap- 



proximately aw 2/3 is valid. The experimental masses 
are reproduced wih an accuracy better than 0.1%. 

Extending the standard charmonium spectrum, three 
new particles have been predicted, two of them associated 
with S°(2455), a charmed baryon and F(4260), observed 
recently. The third particle, labeled X<50> according to 
the proposed SO" (3) level scheme, with a predicted mass 
of 4965 ± 10[MeV], has not been experimentally verified 
yet. 

Summarizing the results of our considerations, the pro- 
posed fractional non relativistic Schrodinger type wave 
equation is a powerful alternative for a discussion of char- 
monium properties and extends our knowledge beyond 
the standard achieved with phenomenological models. 

Therefore fractional wave equations may play an im- 
portant role in our understanding of particles with quark- 
like properties, e.g. confinement. 
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